(N 

O 

(N 



00 



Oh 



> 



Mirror symmetry breaking with limited enantioselective autocatalysis and 
temperature gradients: a stability survey 

Celia Blanco, ^'B Josep M. Rib6,^'^'[l| Joaquim Crusats,^'^'[|| Zoubir 
El-Hachenii,2'3'[l Albert Moyano,^'^ and David Hochbergi '0 

'Centra de Astrobiologia (CSIC-INTA), Carretera Ajalvir Kilometro 4, 28850 Torrejon de Ardoz, Madrid, Spain 
Department of Organic Chemistry, University of Barcelona, 
c. Marti i Pranques 1, 08028 Barcelona, Catalonia, Spain. 
"^Institute of Cosmos Science, University of Barcelona (lEEC-UB), 



c. Marti i Franques 1, 08028 Barcelona, Catalonia, Spain. 

, We analyze limited enantioselective (LES) autocatalysis in a temperature gradient and with in- 

ternal flow/recycling of hot and cold material. Microreversibility forbids broken mirror symmetry 
O ' for LES in the presence of a temperature gradient alone. This symmetry can be broken however 

O , when the auto-catalysis and limited enantioselective catalysis are each localized within the regions 

. of low and high temperature, respectively. This scheme has been recently proposed as a plausible 

model for spontaneous emergence of chirality in abyssal hydrothermal vents. Regions in chemical 
parameter space are mapped out in which the racemic state is unstable and bifurcates to chiral 
solutions. 



INTRODUCTION 



' Recent experimental reports on the deracemization of racemic mixtures of crystals and on the crystallization from 
^ ' boiling solutions [H-Q are striking examples of novel scenarios for spontaneous mirror symmetry breaking (SMSB) for 
<■ compounds for which the homochiral interactions are favored over the heterochiral ones. In other words, these are 
reactions that cannot be explained by Frank-like mechanisms, in which the heterochiral interaction is the favored one. 



^ . Despite some controversy about the actual mechanisms responsible for the SMSB in these situations, the experimental 



reports all coincide in that the final state is stationary: mechano-stationary in the case of wet grinding of racemic 
mixtures of crystals IsHsl and the presence of temperature gradients in the case of deracemization and crystallization 
in boiling solutions [l|, |2|. The above reports stress the fact that the racemic conglomerate crystal mixtures are 
deracemized in experimental conditions where chemical equilibrium is not possible, i.e., they require specific energy 
input to only some of the species of the system (as in the crystal grinding experiments) or else a non-uniform 
temperature distribution. 

The point of departure for the present paper is: if limited enantioselectivity in experimental conditions of closed 
systems with a uniform distribution of temperature and energy inexorably yields a final racemic state, Q then can 
SMSB occur in scenarios of non-uniform temperature distributions? As it turns out, this is in fact possible in the 
case that, in addition to a non-uniform temperature distribution, the systems possess a compartmentalization of the 
enantioselective and the non-enantioselective autocatalyses. 
psj ' The emergence of chirality in enantioselective autocatalysis for compounds which do not follow Frank-like schemes 
T-H I is investigated here for the limited enantioselectivity (LES) model composed of coupled enantioselective and non- 
^ . enantioselective autocatalyses. The basic model [Tj is defined by the following chemical transformations. Production 
of chiral compounds L,D from an achiral substrate A: 

^ ■ 

?H ■ fcl fcl 

_C^_. A . - L, A . - D. (1) 

k-i fc_i 
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Autocatalytic production: 



Limited enantioselectivity: 



k2 k2 
A + L T L + L, A + Dt ^D + D. (2) 



A + L . - L + D, A + Dt ^D + L. (3) 



In contrast to the Frank model, LES is able to account for two important facts: namely, (i) the enantioselectivity of 
any chiral catalyst is limited because of the third reaction Eq.(p]), and (ii) the kinetic link between mirror conjugate 
processes arises from the reversibility of the catalytic stage l'^. The inverse reaction of the non-enantioselective 
autocatalysis (reaction ([3])) substitutes for the mutual inhibition reaction in the Frank model or formation of the 
heterodimer (L + D — > P). Earlier reports had claimed spontaneous mirror symmetry breaking (SMSB) in LES, but 
this cannot occur in either open or closed systems with a uniform temperature distribution. The obstacle comes from 
microreversibility, where K{T) is the temperature dependent equilibrium constant: 

A=i^(r), (l<^<3). (4) 

The condition for the instability of the racemic state is that 

where g = and w = ^ J^. From Eq. (O, we must have 1 - it; > so that 1 > w. But from Eq. dU, || = 
which is incompatible with 1 > w and g < I. This is the situation for open systems. For closed systems, it can be 
shown Q that the racemic state is unstable provided that 

o<9< g.tr' < 5cTr = < 1, (6) 

where the critical parameter in a closed system is always bounded above by the corresponding one for open systems, 
and approaches the latter from below in the limit of large total concentrations C, that is: 9^1°^^'^ gcrtt^- ^1- ® 
is also incompatible with Eq. Q. So the racemic state is always asymptotically stable in this scheme for both open 
and closed systems held at a uniform temperature. Therefore there can be no asymptotically stable chiral outcome 
in this model. 

Nor can the LES model lead to SMSB in closed systems even with a stationary non-uniform temperature distribution. 
On the other hand, numerical simulations of chemical kinetics in a two- compartment model 9] demonstrate that 
SMSB may occur if both catalytic reactions ([2]) and ([3]) are spatially separated at different temperatures in different 
compartments but coupled under the action of a continous internal flow. In such conditions the system can evolve, for 
certain reaction and system parameters, towards a chiral stationary state, i.e., the system is able to reach a bifurcation 
point leading to SMSB. 

This is an appealing result since numerical simulations using reasonable chemical parameters suggest that an 
adequate scenario for such an LES-based SMSB would be that of abyssal hydrothermal vents, by virtue of the typical 
temperature gradients found there and the role of inorganic solids mediating chemical reactions in an enzyme-like 
role. We therefore proposed [1] that a natural prebiotic scenario for such emergence of chirality is that of abyssal 
hydrothermal vents and volcanic plumes [l0l - IT3| which do have the adequate temperature gradients and contain solids, 
as for example clays, which have been proposed by several authors |14l4l8| as catalysts in the prebiotic synthesis of 
organic compounds. 

In view of the above, this paper deals with an analytic/numerical study of the conditions leading to the instability 
of the ideal racemic composition for the LES model with compartmentalized catalysis ([2]) and ([3]) in regions held 
at different temperatures. The two-compartment model (Sec|TT| is already sufficiently involved as to make deriving 
general analytic stability results a near impossible task. We thus focus our efforts on analyzing properties of the 
racemic fixed point; the analytic conditions for its linear stability can be set up and then tested in numerical domains. 
The direct study of the stationary chiral solutions is substantially more complicated and algebraically unwieldy. 
Nevertheless, the characterization of a racemic state as unstable necessarily implies that the system evolves to a 
state of non-racemic composition. Such a nonracemic state could be stable, chaotic or oscillatory. However, the 



3 



many numerical SMSB tests performed for reasonable reaction parameters and under the thermodynamic constraints 
imposed by the principle of microreversibility and the temperature dependencies of the reaction rate constants (see 
Sec lIIIj) have in all cases led to a stable chiral state. These stable states are scalemic mixtures of enantiomeric excesses 
(ee) whose values depend on the parameters of the phase diagram. We first review the impossibility of SMSB in LES 
with a stable temperature gradient. We then prove that SMSB in such a system is possible when the two catalyses ([2]) 
and ([3|) are compartmentalized (localized) in different temperature regions connected by an internal flow of material. 



II. 



LES WITH A TEMPERATURE GRADIENT 



Consider the LES scheme in a temperature gradient: this permits the reaction rates to vary spatially in the system, 
from one place to another, and might provide a way to achieve mirror symmetry breaking. The reverse reaction of 
([3]) in one region (compartment) could be faster that the reverse of ^ in the other. And then mixing could bring 
the hot and cold material into contact. We model this by a closed two-compartment system with volumes V,V*; 
each compartment held at a uniform temperature, T* > T, and internally coupled by a constant internal flow or 
recycling. The concentrations, and reaction rates for the second compartment are labeled by an asterisk, see Figure [1] 
We extend the basic LES scheme Eqs. (|ll2l3l) . by including the terms corresponding to the internal flow. The internal 



(l)+(2)+(3) 



[A*], [D*], [L*] 



(l)+(2)+(3) 



V, T 



[A], [D], [L] 



FIG. 1: Limited enantioselectivity (LES) in two compartments with vohimns V and V* each held at different temperatures 
T* > T and interconnected by an internal flow a of material. See Ref. [9j]. 

"flow" parameter a has units of volume/time; V and V* denote the volumes of each compartment. The corresponding 
rate equations for the two-compartment system have been derived in Q. 

A. New variables 

For setting up a stability analysis it is convenient to employ the sums and differences of the concentrations x — 
[L] + [D], y ^ [L] - [D], X* ~ [L*] + [D*] and y* — [L*] — [D*]. In terms of these, the rate equations take the following 
form: 

M = ^2kM]-iHA] + h[A]-k-i)x + ^{y' + x') + '^ix'-y') + ^{[A*]-[A]), (7) 



^ = 2k,[A] + iHA] - fc-i)x + h[A]x - fc-2 ^ - fe-3 ^ + f (X* - X), 



(8) 



^ - {k2[A]~k^i)y-k^2Xy-k3[A]y + ^{y* -y), 



(9) 



4a . -.:,.,-,,..H.,m-.:..-..w'..-^).^^(x--.v^fl.,-m,, m 
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-j^ = 2kt[A*] + {k;[A*] ~ k*_,)x* + k*,[A*]x* - kU ^^ y ' - fcl3 ^^ ' +y^{x~ X*), (H) 

^ = ik;[A*] - fc* - kUx*y* - k*,[A*]y* + ^{y - y*). (12) 
They satisfy the constant mass constraint: 

V{[A]+x) + V*{[A*]+x*) = C, (13) 
where C is the total conserved mass in the complete two-compartment system. 

III. CONSTRAINTS FROM ARRHENIUS-EYRING 

Certain thermodynamic relationships hold among the reaction rates in both compartments. These will be used 
to prove that SMSB is also impossible for the scheme presented in Sec |lT] in a background temperature gradient. 
Following this demonstration, we then introduce the variant of LES that can and does lead to SMSB. 

From Arrhenius-Eyring, the forward (and reverse) reaction rates for reaction i at temperatures T* and y < y* are 

M 



( )e UT , = (— -)e nT ^ 



(14) 
(15) 



here 

AG,{T)^AH,~TAS, (16) 

denotes the difference in free energy between the activated state (transition state) and the reactants, while 

AG_,(T) = Ai/_, -TA5_„ (17) 

is the free energy difference between the activated state (transition state) and the products; H and S denote the 
enthalpy and entropy, respectively. From the above we can obtain a relation between the forward reaction rates i in 
each compartment: 



^.-(^)exp(-^(^-i))fc,. (18) 



Clearly, once the values of the ki are chosen for the reference compartment at T, we are not free to independently 
choose the reaction rates k* at the higher temperature T*. 

The fundamental microreversibility condition Eq.Q together with Arrhenius-Eyring implies 



ki e 



_ g(AG_.-AG.)/flT _ 



K{T), (19) 



e RT 

^ AG_, - AGi = AAG, (20) 

that is, the difference of the free energy differences AAG must be independent of i, that is, independent of the specific 
ith reaction. This implies that the individual double differences in enthalpy and in entropy must also be independent 
of reaction i, so we must also have 

AAH = {AH_,- AH,), (21) 
AAS'= (AS*., - AS*,). (22) 
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This also gives us an expression for calculating K: 

K{T)^e^e-^. (23) 
The inverse reaction rates are obtained through the constraint Eq.(|4]) as follows: 

We also note that if the constraint Eq. is satisfied at one specific temperature T, then it will automatically hold 
at all others, that is 

K = K{T*) = e^e-^, (1 < z < 3). (25) 

The ratio of the equilibrium constants is given by 

K{T*) /AAi/ , 1 1 ^^ 

The algebraic intricacy of the model in Sec |TT] is already such that we are unable to obtain useful and manageable 
analytic closed form expressions for the conditions leading to the instability of the racemic solution. The situation is 
even worse for obtaining analytic information regarding the possible stationary chiral solutions. We appeal instead to 
chemically inspired conjectures that can be tested numerically for coherence and compatibility with microreversibility. 

First, it is clear that in view of the gradient T < T*, the putative condition, which could conceivably lead to 
symmetry breaking in the limit of small values of a, 

fc_2 < & h > h at T 



and < ^-3 ^ ^2 > ^3 ^* 



(27) 



is incompatible with the constraints in Eqs. (|23I25|) . This condition (l?7l) is inspired by the observation that for a 0, 
the two compartments are practically isolated from each other and can be treated as approximately independent. 
These are thus the necessary conditions for obtaining an unstable racemic solution in each compartment (see Sec HI . 
But they are incompatible with microreversibility. 

Secondly, the analysis in Sec U suggests that SMSB might occur when the inverse reaction of ([3]) in one region is 
faster than the inverse reaction of ([2]) in the other region. Taking microreversibility into account, the only way this 
might be achieved is, for example, by arranging for 



fc_3 < fc_2 & ^2 > fca at T 
and k*_2 < fc* 3 & ^2 < ^3 at T* 



(28) 



But this is forbidden by virtue of Eq. (1261) . which is satisfied by the ratio of the equilibrium constants. So a 
temperature gradient and internal flow are by themselves, not enough to produce a bifurcation. Actually, no spatially 
varying temperature profile is sufficient, as can be seen by partitioning the closed system into a number of sufficiently 
small regions within which the local temperature is approximately uniform. 



IV. TEMPERATURE GRADIENT AND IMMOBILIZED CATALYSTS 



Our working hypothesis is that a necessary but not sufficient condition for the instability of the racemic solution is 
k^3 > k-2 and /c2 > ^3 in the presence of immobilized catalysts that ensure that kj_2 = in one region and k±^ — 
in the other Q. Other conditions (total system concentration C, the flow rate a, compartment volumes V, V* , etc.) 
also come into play for determining the overall instability, in a highly nontrivial and nonlinear fashion. 
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A ■ 



( (-2fci - [fc2 + HX - f - (-(fe + fc3)A + (fc-2 + fc-3)x + - 

(2A:i + [fc2 + fcalx) ((fc2 + fcs)^ - fc-i - (fc-2 + fc-3)x - f ) f 

+ + fci]x*) - ^(2fc! + [fci + ki\x') m + - ^(.A + x) - r ) 

-fell - (fcl2 + fcl3)X* - T?^ 

-(2fci* + [fc2* + fc3]r)) / 



(31) 



A. Linear stability analysis of the stationary racemic fixed point 

The equations (|9ll2p for ^ = = are identically satisfied for the stationary solution y = y* — 0. We therefore 
carry out a stability analysis of this racemic fixed point and determine whether the racemic solution (racemic in 
both compartments) is asymptotically stable or unstable. This will depend on the internal flow parameter a that 
characterizes the cycling of hot to cold material between the two compartments. Clearly, if we set a = 0, we merely 
recover two isolated copies of LES in independent closed compartments, each one at a constant temperature, and 
there can be no mirror symmetry breaking in this situation; the considerations of Sec |T] apply. 

An algebraic advantage of studying the racemic fixed point y = y* — is that the five independent concentration 
fluctuations decouple into two sets of three and two, respectively. This situation is reflected in the structure of the 
Jacobian J which then reduces to a block-diagonal form with a 3 x 3 sub-block corresponding to the fluctuations 
{6A,5x,Sx*) and a 2 x 2 sub-block corresponding to {Sy,Sy*) thus: 

/^3X3 ^3X2 \ /^3X3 q \ 

■J -[^ ^2X3 52x2 J ^ Q 52x2 J ■ (2J) 

The temporal evolution of the linearized concentration fluctuations about the racemic fixed point y = y* = of the 
kinetic equations is given by 



where the 3x3 array A is given by Eq. (|3ip and the 2x2 array B is given by Eq. ([32)) , (see Appendix 
The Jacobian matrix J (1^ must be evaluated on non- negative stationary solutions ^>0,x>0,x*>0 corre- 
sponding to y = y* = 0. Then, to assess the stability of the solution, the five eigenvalues Xi,i — 1,2, ...,5 of the 
Jacobian must be calculated. If any one of these five eigenvalues is positive (or their real part, if complex) then the 
solution is unstable. This means that the system will evolve to a chiral state, so mirror symmetry will be broken. Only 
if all five of the eigenvalues are negative (real part) can we claim that the solution is stable. Deriving manageable and 
useful closed form expressions for the eigenvalues Xi of Eqs. (I31l32p is practically impossible, due to the fact that the 
racemic fixed point solutions A, x, x* l^ad to unwieldy expressions (as solutions of coupled quartic equations). On the 
other hand, direct numeric calculation of the fixed point solutions and their associated eigenvalues is amenable and 
provides a wealth of information about the dynamic stability of the underlying model, as functions of the chemical 
rates and the system parameters. We will therefore map out regions of stability/instability in parameter space. We 
carry this out assuming immobilized catalyst from the start, setting fcji ^-2 to zero in one compartment, and ^3, fc_3 
in the other. Variations in the remaining rate constants are carried out satisfying the constraints in Eqs. (|24l25p . 

We apply a second, independent, stability test which does not require calculating the eigenvalues: namely, the 
Routh-Hurwitz (RH) criteria. We derive explicit expressions whose algebraic signs indicate whether the racemic fixed 
point is stable or unstable. The canonical form of the characteristic polynomial for the complete 5x5 jacobian Eq. 
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(IMl) is 

P{X) = + aiX^ + OiA^ + asX^ + + as = 0. (33) 

Then (see Appendix B of [20]) there are conditions on the coefficients Oi, i = 1,2,..., 5 such that the zeros of P{X) 

have KA < 0. The necessary and sufficient conditions for this to hold are the Routh-Hurwitz conditions. One such 
form, together with 

og = - det (B) det (A) > 0, (34) 

is that 

Di = ai = -[tr(S) + tr(A)] > 0, & (35) 

D2 = {aia2 - as) > 0, & (36) 

D3 = a3D2 + ai(a5 - aia4) > 0, & (37) 

Di = a^Ds - a^laiial - a^) - {a3a2 ~ 05)} > 0. (38) 

The expressions a, can be read off directly from comparing the polynomial in Eg. ([55]) to — P calculated in Eq. (|B5I) . 

If any of the above conditions Eqs. (1341138^ does not hold, then the racemic fixed point solution is unstable. As 
before, this means the system will evolve to a chiral state. This test can be compared with direct numerical calculation 
of the five eigenvalues Xi i = 1,2,3,4,5 (the roots of the characteristic polynomial). We find complete agreement 
between the two methods (eigenvalues, RH criteria) employed. We emphasize that an instability in the racemic fixed 
point implies the onset of a bifurcation to a non-racemic, hence chiral, solution. 

B. Domains of instability 

We initiate the procedure outlined above by specifying the forward/reverse reaction rates (temperature differences 
are treated implicitly), the internal flow rate, the compartment volumes and the conserved total system mass: 

{k±,iT),kUT*)},a,V,V*,C, (39) 

the individual rates of course satisfying microreversibility at the respective temperatures T, T*. We then solve for the 
complete racemic fixed point solution, retaining only those solutions that are non-negative: 

We next evaluate the Jacobian matrix over this fixed point solution: 

/ ^3x3 g \ 

"^=^ ( £2X2 j \{A,x,x'}>o- (41) 



As a final step we evaluate the five RH conditions Eqs. (|34l - 1381) : are they all true or not? If not, then we immediately 
know that the racemic fixed point is unstable for the parameter choice made in Eq. p9p . Hence any fluctuation about 
the idealized racemic composition will grow and drive the system to a chiral final state. In parallel, we also evaluate 
numerically the flve roots of the characteristic polynomial Xi i — 1, 2, 3, 4, 5 and verify agreement between RH criteria 
and the eigenvalues. We use the RH criteria to map out the regions of linear instability for the racemic solution. 

We start with the reaction rates ki = 10"^, = 10"l^ fc* = 10"^ fcl^ = 10"", fc2 = 10^, fc_2 = 10"*, k* = 0, 
fcl2 — 0, ks = 0, fc_3 = 0, A;| = 10^, fc* 3 = 10"^, the volume and flux parameters V = 10, V* = 10, a = 0.1, and the 
initial concentrations [^]o = 10-^ [A*]o = 10'^, [L]o = 10"", [D]o = 10"", [L*]o = 10-" + 10-2i, [D*]o = 10"" as 
employed in Ref @. This point in parameter space was shown to lead to SMSB and subsequent chiral ampliflcation 
from direct numerical integration of the differential rate equations. These specific rate values and system parameters 
were obtained after performing a set of trial and error numerical simulations obeying the condition fclg > fc_2 and 
^2 > ^3 which corresponds approximately to thermodynamically unattainable conditions for systems with a uniform 
temperature and lacking compartments. Once obtained, we can exploit the stability analysis to map out and amplify 
the full domain of instability of the racemic solution. 

We consider how the stability of the racemic fixed point responds to variations of selected pairs of variables about 
this point. In Fig [5] we display the regions of stability/instability as a function of (a) the two compartment volumes 
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V and V* , (b) the flow rate a versus equal compartment volumes V = V* , (c) the total concentration C versus flow 

k* 

rate a, and (d) a = versus fc_2. The graphs in (a) and (b) merely tell us that the ratio of the compartment 
volumes cannot be either arbitrarily large nor small, and that for equal compartment volumes, there is a minimum 
flow rate below which no chiral state can be obtained. Perhaps more surprising are the trends indicated in Fig [2] (c) 
and (d): for a given flow rate there is a critical system concentration above which mirror symmetry cannot be broken 
and in (d), for a given value of /c_2 a bounded region in fclg in which the racemic state is unstable, and this region 
narrows down and pinches off for a critical value of fc_2- We can here appreciate that the conjecture A;^3 > A;_2 is 
necessary, but certainly not sufficient to lead to an instability in the racemic state. We therefore delve further into 
this nonlinear relationship by considering how the instability varies with fclg, k-2 and with the flow rate a. In Fig [3] 
we display a sequence of plots showing how the region of instability grows in area as we scale up the flow rate. For 
a given flow rate, there is a maximum value of k-2 beyond which no value of /c* 3 will lead to an instability, and this 
maximum value scales with a. Below this maximum value, there are always upper and lower bounds on k'^^ between 
which the system is unstable. These bounds also scale with a. In Fig [3] we keep both /clg and fc_2 fixed and vary a: 
increased flow enlarges the region of instability. This can also be appreciated in the three-dimensional plot in Fig 0] 
exhibiting how the region of instability expands in cross sectional area as we scale up the internal flow rate a. 

In Fig [5] we display a sequence of plots showing how the region of instability decreases in area as we scale up the 
total concentration C, over three orders of magnitude. As C is scaled up, the rates k*_^,k-2 decrease in magnitude 
in such a way as to preserve the shape of the region of instability. As before, for a given total concentration, there is 
a maximum value of k-2 beyond which no value of k*_^ will lead to an instability, and this maximum value inversely 
scales with C . Below this maximum value, there are always upper and lower bounds on fclg between which the system 
is unstable. These bounds also inversely scale with C. A three-dimensional plot in Fig IH] indicates the cross-sectional 
area of the domain of instability shrinks as we scale up the total system concentration C Dilute concentrations are 
more favorable for SMSB 



V. SYMMETRY BREAKING AND ENTROPY PRODUCTION 



We illustrate the above considerations by way of two examples. The inherent chiral fluctuations about the ideal 
racemic composition can be modeled by starting with an initial ee below the statistical deviation. In the flrst 
simulation. Figured we begin with [A]o = [A*]o = 1 x IQ~^M, [L]q = [D]q = 1 x lO'^M, [L*]o = 1 x 10"^ + 
1 X lO^^^M, and [-D*]o = 1 x lO^^A/. In the second. Figure [51 we keep the same rates and system parameters 
but start off with different initial concentrations: [A]o = [A*]o = 1 x 10""M, [L]o = [D]^ = 5 x [L*]o = 

5 X lO-'^ + 1 X and [D*]o = 5 x IQ-Hi. 

We calculate and superimpose the entropy production[2l| a on the concentration curves. The entropy production 
is a measure of the dissipation of the out-of-equilibrium process under study. For the initial conditions leading to 
Figure [3 the entropy production suffers a peak at the onset of the induction period, reflecting the fact that the 
catalysis of enantiomers Eq. ^ , and subsequent depletion of achiral substrate, is the most important contribution. 
The production then falls to a small but constant nonzero value and remains nonzero as long as the system is kept 
out of equilibrium. In the case of Figure |8l the entropy production starts off from a large value then drops (at the 
same time scale as the peak in Figure [T]) and subsequently drops down to its minimum nonzero value at the symmetry 
breaking bifurcation. It again remains small but nonzero as long as the system is kept out of equilibrium. 



VI. CONCLUDING REMARKS 



Compartmentalization, together with different temperatures for the enantioselective and the non-enantioselective 
autocatalyses, was shown to be a necessary condition for SMSB. The temperature gradient and internal flow or 
recycling of hot and cold material is the required driving force keeping the system far from equilibrium. On the basis 
of the stability analysis of the racemic flnal state, the results presented here demonstrate the existence of SMSB as 
a closed region (see Figs. 3, 5 and 6) in the phase representation of the reaction parameters. Explicit examples of 
SMSB are displayed in Figures [7] and [8] showing also the time dependence of the associated production of entropy for 
starting from two different initial concentrations. 

The considerations given to the microreversibility condition in Sees U |TT] and IIIII serve to underscore the fact 
that limited enantioselectivity-of and by itself- cannot lead to spontaneous mirror symmetry breaking, neither for a 
uniform temperature (i.e. at an experimental condition required to attain chemical equilibrium), not even for two 
temperatures (i.e. an experimental condition that excludes chemical equilibrium). By contrast, mirror symmetry 
breaking may be possible when in addition to a non-uniform temperature distribution the reactions (ii) and (iii) occur 



9 



in two distinct compartments held at different temperatures. The microreversibihty arguments [8|,|22], used as a proof 
of the correctness of the model, are necessary because we assume a scenario where chemical kinetics can be applied. 

We briefly point out what modifications would have to be made in order to consider spatial temperature gradients. 
First, we need to input the temperature profile T{x,y,z) as a known field. This converts the reaction rates hi into 
spatially dependent functions via the Arrhenius relation. The kinetic rate equations must be replaced by partial 
differential reaction-diffusion-advection equations for the spatial and temporal dependence of the concentrations. The 
background convective flow of hot/cold material could be modeled by a steady hydrodynamical flow- field compatible 
with the compartment boundary conditions. In short, the technical complications would be considerable, but not 
insurmountable. 

The significance of the above results with respect to applied absolute synthesis is somewhat limited due to the 
temperature difference required for achieving the essential inequality k*i^ > fc_2. The point is, a sufficiently large 
difference would require high pressure conditions as well to maintain the media in a liquid state. By contrast, the 
results can be of significance in scenarios of prebiotic chemistry. In this context, the required experimental conditions 
of high temperature gradients and compartmentalization agree with those found in deep ocean hydrothermal vents, 
as we have previously reported Q. Furthermore, the more relevant feature of this LES scenario is its ability to 
exhibit SMSB at very low reactant concentrations. In this respect, the detailed description given here for the effect 
on the SMSB dependence of the minimal internal flow rate -with respect to the compartment volumes at different 
temperatures- and on the reaction rate inequality fclg > fc_2, shows that the variability range for SMSB of these 
parameters increases when the total system concentration C decreases, i.e. the probabilities for SMSB to occur 
increase in a prebiotic scenario. Notice that this overcomes the more important difficulty for a reasonable SMSB in a 
prebiotic scenario, where a specific prebiotic organic compound could be present, but only at very low concentrations. 
This, in spite of the presence of an important fraction of organic compounds, because of the diversity of different 
organic compounds in prebiotic scenarios psj . While LES is ruled out as a scheme for SMSB in experimental 
conditions where thermodynamic equilibrium can be achieved, the results presented here open up the study for SMSB 
in scenarios of non-uniform temperature distributions, specific energy inputs to some species of the system, and 
compartmentalization, i.e. in conditions similar to those of living systems. 
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The linearized fluctuation equations that follow from the kinetic equations Eqs. ((71 IT2]) in Sec. Ill Al are as follows. 
The overbar denotes a stationary solution of Eqs. ([Tl fT^ . 
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Appendix A: Fluctuation equations 



SA = ( - 2fci - [fc2 + h]x 



V 



a 



V* 



a 



)SA+(^- (fc2 + h)A + (fc_2 + fc-3)x + fc-i - 



+ (fc_2 - k-3)ySy 




(Al) 



(2fci + [k2 + k3]x}SA 



+ 



( 



(fc2 + k3)A - fc_i - (fc_2 + fc-3)x -^)Sx + {k-3 - k^2)y Sy + ^6x* (A2) 



= {k2 - k3)y6A - fc_2y^X + 



k-i - k^2X + (^2 - ^'3)^ 



(A3) 



Sx* 



= -^^2kl + [k* + ki]r)SA + - ^^2kl + [k* + fcgir)^^ 
+ (fe* + - -fx) - r) k*_, {ku + ku)r 

+ (fc*3 - k*_2)y* 8y* 



{2kl + [kl + kl]r)yx* 



(A4) 



10 



+ {ik;-k;){^~^iA + x)-r)^kl,^kl,r-^)Sy*. (A5) 

Specializing to the racemic fixed point y — y* = leads to the decoupling of the first three fluctuations from the 
latter two: 

SA = ( - 2fci - [k2 + k3]x - - ^)6A + ( - (fc2 + k3)A + (fc_2 + k.3)x + k-, -^)Sx- ^Sx* (A6) 

S'x = (2fci + [k2 + ksWSA + ((fc2 + ks)A ~~ - (fc_2 + k^3)x -^)^X+ ^Sx* (A7) 

S'y = + ( - fc-i - fc-2X + (fc2 - ks)A -^)Sy + ^Sy* (A8) 

^x* = -^m + [k; + k;]x*)SA+[^-^{2ki + [k; + k;]r)sx 



{[k; + k;]{^ ~ ^{A +x)- xi k*_, - {ku + ku)r - - m + + ki]r)yx* 



(A9) 



as reflected by the specific matrix entries in Eqs. pil32p . 

Appendix B: Characteristic polynomial 



The characteristic polynomial associated with the Jacobian matrix Eq. (1291) evaluated at the racemic fixed point is 

P(A) = det|J- A/| =det|A- A/|det|B- A/|, 

= P(3'(A)P(2)(A)^ (Bl) 

where / is the identity matrix and the quadratic and cubic polynomials are given by 

p(2)(A) = A2-tr(P)A + det(P), (B2) 
P(=^)(A) = -A^ + tr(A)A2 +G(A)A + dct(A), (B3) 

respectively, where 

G[A) = -aiia22 - anasa - a22a33 + a32a23 + ai2a2i + ai^a^i- (B4) 

Thus, inserting Eqs. (|B2IB3|) into Eq. ((BT|) . we obtain the fifth order polynomial 

P(A) = -A-V [tr(P) +tr(^)]A* + {G(A) -det(P) -tr(P)tr(^)}A^ 

+ {tr(A) det [B] - G(A)tr(P) + dct {A)}\^ + {G{A) det (P) - det (v4)tr(P)}A 

+ det (P) det (^) . (B5) 
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FIG. 2: White zones indicate where the racemic state is unstable to perturbations and so bifurcates to a chiral state 
(color/grayscale, where racemic state is stable) as a function of the indicated selected pairs of variables: a) The cold and 
hot compartment volumes V and V* , b) the flow parameter a versus the volume of each compartment, holding V = V* fixed, 
c) the flow parameter a and the total system concentration C, and d) the reaction rate k-2 and a, where a = fe* 3/^-2. Except 
for the specific pair that is varied in a),b),c), and d), the remainder of values are held fixed at fci = 10"^, fc-i = 10"^^, 
kl = 10"^, fell = 10"", k2 = 10^ k-2 = 10"*, k2 = 0, fc*2 = 0, = 0, k-3 = 0, k^ ^ 10\ fclg = 10"^, the volume and 
flux parameters V = 10, V = 10, a = 0.1. The initial concentrations [A]o = lO^'^, [A*]o = 10"'', [L]o = 10"", [D]o = 10"", 
[L*]o = 10"" + 10"2i, [D*]o = 10"". 
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FIG. 3: White zones indicate where the racemic state is unstable to perturbations (color/grayscale, where its stable) for 
different values of k-2 and fc* 3, and for different values of the internal flow parameter a. The area of the region of instability 
scales up with the flow rate. The remainder of parameters and initial concentrations as in Fig. [2] 
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FIG. 4: Effect of varying the internal flow rate. Three dimensional figure represents all points for which the racemic state is 
unstable to perturbations. Increased flow enlarges the allowed region of instability. The racemic state is stable in the empty 
(white) domain. The remainder of parameters and initial concentrations as in Fig. [5] 
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FIG. 5: Regions where the racemic state bifurcates to chiral states (in white) for different values of fc_2 and fc*3, and for 
different values of the total concentration C. Increasing C decreases the ranges in both k-2 and fc* 3 for which the racemic 
state is linearly unstable. The remainder of parameters and initial concentrations as in Fig. [2l 
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FIG. 6: Effect of varying the total system concentration C. The soUd three-dimensional figure (color/gray) indicates where 
the racemic state is unstable to perturbations. Increasing the system concentration C shrinks the cross sectional domain 
k-2,k*_3 of instability. The racemic state is linearly stable in the empty (white) region. The remainder of parameters and 
initial concentrations are as listed in Fig. [21 
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FIG. 7: Symmetry breaking bifurcation in the two compartment model, see Figure [T] For the rates and system parameters in 
given in Sec. IIVBI The entropy production a peaks at the onset of the induction period, well before the symmetry breaking 
bifurcation. Note cr > for all times. 
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FIG. 8: Symmetry breaking bifurcation in the two compartment model, same rates and system parameters as in Figure [71 See 
text for the initial concentrations. In this case the entropy production a (scaled up by a factor of lO'^) decreases monotonically 
and drops to a minimum (but nonzero) value at the symmetry breaking bifurcation, and remaining positive for all subsequent 
times. 



